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Abstract 

In the Next-to-Minimal SuperSymmetric Model (NMSSM) the lightest CP-odd Higgs bosons (oi) 
can be very light. As a consequence, in addition to the standard charged Higgs boson (h^) decays 

£N( considered in the MSSM for a light charged Higgs (m h ± < m t ), the branching fraction for — > a{W 
can be dominant. We investigate how this signal can be searched for in ti production at the Large 
Hadron Collider (LHC) in the case that (m ai > 2me) with the a\ giving rise to a single 66— jet and 
discuss to what extent the LHC experiments are able to discover such a scenario with an integrated 
luminosity ~ 20 fb _1 . We also discuss the implications of the possible Higgs- signal observed at the 

<C LHC. 
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I. INTRODUCTION 



With the successful start-up of the LHC and the intriguing results from the first year of 
data-taking at the center of mass energy of 7 TeV and a integrated luminosity close to w 5 fb _1 
for each of the ATLAS and CMS experiments, the ongoing run in 2012 is set to be a milestone 
in particle physics. The possible signal for a Higgs boson around 125 GeV may or may not 
be confirmed and the search for physics Beyond the Standard Model (BSM) will continue. 
Contrary to the neutral Higgs boson, which if discovered may need to be analyzed in detail 
regarding its branching fraction into various channels in order to determine whether it is the 
Standard Model (SM) Higgs boson or not, the discovery of a charged Higgs boson would be an 
unmistakable sign of BSM physics. 

The charged Higgs boson arises in theories with more than one Higgs doublet. The prime 
example is the MSSM [TJ [2] which has two Higgs doublets, leading to two CP-even Higgs 
bosons (/i, H) and one CP-odd (A) in the case of CP-conservation and two charged states 
(i? ± ) after electroweak symmetry breaking. In this case the two Higgs doublets are required by 
supersymmetry with one of the giving masses to the up-type fermions and one to the down-type 
ones. For a complete introduction to the Higgs sector in the MSSM we refer to [3]. 

The main reason for introducing supersymmetry (for a general introduction to supersym- 
metric theories we refer to [4J) is to solve the so called hierarchy problem, i.e. why the scale of 
the electroweak interaction is so much smaller than the cut-off of the SM, normally taken to be 
the Planck mass where gravity becomes the dominant force and the SM breaks down. With the 
Higgs boson being a scalar particle, the higher order corrections to its mass are proportional 
to this cut-off and with the SM only being an effective theory this cut-off dependence cannot 
be renormalized away. In a supersymmetric theory this problem is essentially solved by the 
introduction of the fermionic partners of the Higgs fields, which only get logarithmic corrections 
to their masses and thereby avoids fine-tuning. In turn this means that the Higgs boson masses 
are also protected from receiving quadratic corrections as long as supersymmetry is not broken 
or only softly broken. 

In addition to solving the hierarchy problem supersymmetry also offers a candidate for cold 
dark matter [5j[6] in the case that i?-symmetry is preserved, which in turn is introduced to avoid 
terms in the Lagrangian that otherwise would mediate proton decay. In this case the lightest 
supersymmetric particle (LSP) is stable and generically it has the right mass and cross-section 
to constitute the observed dark matter. Supersymmetry also improves the unification of gauge 
forces at an hypothesized grand unification scale although the unification is not exact and it 
does depend on the details of the spectrum of SUSY-particles. 

As is well known, the MSSM by itself is not without problems. Leaving the question of the 
precise mechanism for supersymmetry breaking aside, the MSSM faces the so called //-problem. 
This relates to the magnitude of the dimensionful /x-parameter which couples the two Higgs 
doublets to each other in the superpotential. In order to avoid large cancellations between this 
contribution to the Higgs masses and the soft supersymmetry breaking terms as well as having 
a phenomenologically viable supersymmetric theory (mainly having a large enough chargino 
mass), the magnitude of \± should be of order the electroweak or supersymmetry breaking 
scales. The problem is then that there is no a priori reason for this parameter to have any 
particular value, in principle it could be anything up to the Planck scale, so why is it similar 
to the electroweak or supersymmetry breaking scales? 

In the NMSSM the /x-problem is solved by introducing an additional Higgs singlet into the 
theory. After supersymmetry breaking this field gets a vacuum expectation value (vev) that 
effectively acts as a /i-term. The original /x-term in the superpotential can then be set to zero 
without spoiling the viability of the theory. For a more detailed review of the NMSSM we refer 
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to [HE!. 



The additional Higgs singlet has important consequences for the phenomenology of the Higgs 
sector. In the MSSM, the masses of the heavy Higgs bosons (H^A^H^ are closely related to 
each other as they originate from the same (second) Higgs doublet if viewed in the Higgs basis 
where only one (the first) of the Higgs doublets has a vev. For example, at tree-level the masses 
of the CP-odd and charged Higgs bosons are related by m 2 H ± = m 2 A + m^. As a consequence 
the decay H ± —> AW is typically not open. 

In the NMSSM the additional Higgs singlet means that there is one more CP-even and 
one more CP-odd field with a separate mass scale introduced into the Higgs sector. As a 
consequence, the by now three CP-even and two CP-odd electroweak states will mix into the 
respective mass eigenstates. Thus the mass-relations from the MSSM will be altered. This is 
particularly evident in the CP-odd sector where the lightest state a\ may now be much lighter 
than the charged Higgs boson - even after taking experimental constraints into account as 
discussed below - opening up the possibility for the a\W decay to be dominant. In turn 

this means that the search for charged Higgs bosons, in t-quark decays for example, has to be 
widened also to include this decay channel. 

The decay —> a\W has already been considered to different levels of detail [9l-fTTj in the 
literature and there are constraints from the DELPHI experiment for m ai > 12 GeV [12], as 
well as the CDF experiment for the case a\ —± t + t~ [13]. In this paper we want to focus on 
the region in parameter space where the a\ mass is above the bb threshold but still so close to 
it that the two 6-quarks will fragment into a single 66— jet. The viability of scenarios with light 
ai's have also been considered by [HHT8]. 

Our paper is organized as follows. In the next section we give some basic properties of 
the Higgs sector in the NMSSM that are relevant to our discussion. We then discuss the 
constraints on the parameter space in section 3, including the latest results from LHC. In 
section 4 we illustrate how the signal a{W can be searched for in it-production taking 

into account the appropriate backgrounds. Section 5 contains a discussion of the implications of 
the possible Higgs signal from the ATLAS and CMS experiments and in section 6 we summarize 
and conclude. 



II. BASIC PROPERTIES OF THE NMSSM 



We consider the Z 3 -symmetric version of the NMSSM with the superpotential given by 

V^NMSSM W^MSSM + \SH U • + (l) 

where Wmssm is the superpotential of the MSSM with /i set to zero. The soft supersymmetry 
breaking potential relative to the MSSM is then given by 

y NMSSM = K MSSM + ^2 + . ^ + 1^ + h ^ (2) 

where the part of V^ft SM only depending on the Higgs fields is given by, 

T/ MSSM _ 2 |zt|2| 2 I tt |2 fn\ 

Ksoft^iggs - m H u \ H u\ +m Hd \H d \ . [6) 

In addition SM contains all the dependence on the other soft supersymmetry breaking 
parameters: the gaugino masses Mi,M 2 ,M 3 , the tri-linear couplings a w ,ad,a e , the squark 
masses Mq, M w , M^, and finally the slept on masses M^, M e . In the following we will assume 
minimal flavour violation so that the sfermion mass matrices are diagonal and the tri-linear 
couplings are proportional to the corresponding Yukawa coupling matrices r u = A u y u etc. 
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After electroweak symmetry breaking, and assuming that CP is conserved, the Higgs sector 
will contain three CP-even Higgs bosons (hi, h 2 , ^3), two CP-odd (ai, a 2 ) and one charged (/i ± ), 
where the states are ordered in terms of increasing mass. In the same way as in the MSSM the 
minimization conditions for the Higgs potential allows one to trade the vrtu u , rnn d parameters 
for the doublet vev v ~ 174 GeV and tan/3 = v u /vd- Similarly ms can be expressed in terms 
of the singlet vev v s which in turn gives rise to the effective /x-parameter, \i = Xv s . All in all 
this leaves us with 6 unknown parameters describing the Higgs sector of NMSSM at tree-level: 
tan/5, /i, A, ft, A\, A K . Below we will trade the latter two parameters for the masses m h ± and 
m ai . 

As already alluded to the mass-eigenstates are mixtures of the electroweak eigenstates. 
More specifically, writing S weak = (Re(# u ), Re(H d ), Re(S)) we have h % = S ZJ SJ eak . (In the 
MSSM limit this means that Si 2 = — sin a.) Similarly for the CP-odd states we have = 
AijAj edk with A weak = (Im(cos/3i7 w + sin^iifd), Ini(S'))- Here the mixing matrix is simply 

A = ( COS ^ S ^ n ^ n A ) . Together with the ratio of the two doublet vevs (or equivalently the 
y - sin 9 a cos 9 a J 

rotation angle needed to go to the Higgs basis where only one of the doublets have a vev) tan /3, 
the mixing matrices S and A specify the reduced couplings to fermions and gauge bosons as 
given in table |TJ 



Vertex NMSSM MSSM SM 



hitt 


Sn 


cos a 


1 


sin/3 


sin/3 


hibb 


S12 


sin a 


1 


cos f3 


cos /3 


h 2 tt 


S21 

sin (3 


sin a 
sin /3 


n.a. 


h 2 bb 


S22 


cos a 


n.a. 


cos ft 


cos /3 


a\tt 


cot (3 cos 9 a 


cot /3 


n.a. 


a\bb 


tan/3 cos6U 


tan /3 


n.a. 


h±VV 


sin /3 Sn + cos /3 S12 


sin(/3 — a) 


1 


h 2 VV 


sin/3S 2 i + cos/3S 2 2 


cos(/3 — a) 


n.a. 


a\h\Z 


(cos /3 Sn — sin /3 S12) cos 9a 


cos(/3 — a) 


n.a. 


a\h 2 Z 


(cos /3 S21 — sin /3 S22) cos 9a 


sin(/3 — a) 


n.a. 


hih + W~ 


cos /3 Sn — sin /3 S12 


cos(/3 — a) 


n.a. 


ai h + W~ 


cos 9 a 


1 


n.a. 



TABLE I. Reduced Higgs couplings in the NMSSM compared to the MSSM and the SM (when 
applicable). Note that the reduced couplings to fermions are identical for all three generations, even 
if only the third generation is displayed here. The couplings to /13 can be obtained from the hi ones 
by the replacements Sn S31 and S12 — > S32 wheras the couplings to a 2 can be obtained from the 
ai ones by the replacement cos 9 a sin 9a- 



III. EXPERIMENTAL CONSTRAINTS 



In this section we will explore to what extent the process we are interested in is constrained 
by existing experimental data. Since we are interested in a light (with m^± < m t ) and 
a light a>i there are constraints both from collider experiments as well as low-energy flavour 
experiments. However, before going in to the various constraints we will specify the scenarios 
that we have considered and then come back to the question of experimental constraints. 
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A. Specification of SUSY scenario considered 



In the following we will consider a variant of the well motivated m™ ax -scenario in the MSSM 
similarly to what was done in [T7j . Thus we will consider a universal scale Msusy for the 
sfermion masses at the supersymmetry breaking scale. In other words we assume, as already 
stated, that the sfermion mass matrices (Mq etc) are diagonal, and furthermore we assume that 
all diagonal entries are equal to Msusy which we keep fixed at 1 TeV. In addition we assume that 
the gaugino masses are related as in the constrained MSSM, where supersymmetry breaking is 
assumed to be mediated by gravity, namely M x = 100 GeV, M 2 = 200 GeV, M 3 = 800 GeV. 
Finally we will assume that A t = A^ = A T but contrary to what was done in [17J we will let 
them vary in the range A t £ [—5000, 5000] GeV so that the amount of mixing between the II 
and Ir is unconstrained. 

For the Higgs sector we will let all six parameters vary freely. However, as was done in [T7j, 
we will trade the A K parameter for m ai and A\ parameter for m^± using an iterative procedure 
starting from the tree-level relations: 



M 



2 



A / A k \ \ 2 v 2 sin 2/3 / A t k \ k A 



(5) 



2/i 



where the latter gives the masses of the mass-eigenstates ai, a 2 after diagonalisation. Thus the 
parameters we consider with their respective ranges are: 



tan/3 e [1,60], 
A e [0,0.7], 
k e [-0.7,0.7], 
/x e [125, 1000] GeV, 
m h ± e [80, 170] GeV, 
m ai e [4, 150] GeV. 



The limits for the various parameters has been chosen as follows: For tan /3 , k and A we impose 
perturbativity up to the GUT scale which effectively means that any value out side the above 
regions are bound to fail. (In addition some points inside these regions also fail because of this 
requirement.) The lower limits on \± and rrih± are dictated by experimental constraints. The 
upper limit on ji is not a hard one but follows from the implicit assumption that /i should be 
of order the electroweak scale whereas the upper limit on m^i is given by the condition that 
the decay t — >> bh + should be open. The reason for letting /i vary freely is mainly that this 
decreases the correlations between the masses of the Higgs bosons as will be discussed more 
below. Finally the lower limit on m ai is chosen in order to have a\ — >■ r + r~ open whereas the 
upper limit follows from having m h ± < 170 GeV. 

In order to calculate the resulting models from the inputs we use the package NMSSMTools 
version 3.2.0 [20l I2T] with default settings. Among other thing this means that we impose 
perturbativity of the model up to the GUT scale. Finally, in all scans we generate ~ 1 M 
points (with flat priors in the parameters considered) which fulfill the theoretical constraints 
implemented in NMSSMTools. 
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B. Current experimental constraints 



The most important constraints comes from the direct searches for Higgs bosons for which 
we use the package HiggsBounds version 3.7.0 [221 E3J. In addition there are also constraints 
from direct searches for supersymmetric particles, various flavour constraints and in principle 
also the anomalous magnetic moment of the muon as well as the relic density of dark matter. 
We have not applied the lattr two constraints for the following reasons. To investigate the 
amount of dark matter in the various models one would also need to vary the gaugino masses 
as was done in [24J. Since we keep these fixed we have not applied the dark matter constraints. 
On the same vain we have not applied the constraint from the anomalous magnetic moment 
of the muon, since this depends on the masses of the scalar partners of the muon and the 
neutrinos, apart from requiring /i to be positive. 

When it comes to the flavour constraints the situation is more involved in that various 
constraints have different level of model dependence. On the one hand there are constraints 
from tree-level mediated processes such as B u —> r + v T) which only depend on the Higgs sector 
and on the other hand there a constraints from loop-mediated processes such as B s — >> and 
b —> 57 which depend on details of the supersymmetric sector of the model. In the following we 
will limit ourselves to applying the most severe constraints from B u r Jr v T) which limits the 
available parameter space in [m h ± , tan /3] and B s —> which puts limits on [m ai , tan /3] . The 

latter constraint is especially important since we will consider light a± which is very constrained 
by the data from LHCb and CMS [25], [26]. Finally we apply the direct constraints from searches 
for supersymmetric particles. For all constraints except the Higgs bosons we use NMSSMTools 
version 3.2.0. 

The results from the scan are displayed in Fig. [T] with black points being viable models 
from a theory point of view but excluded by the direct searches for Higgs bosons and coloured 
points being allowed by the same constraint. Of primary interest are the allowed regions in 
[m ai ,m/j±]. As is clear from the figure there is a distinct region of points with m ai ~ 2rriB ~ 
10.6 GeV allowed by all constraints (indicated by green colour) for essentially any value of 
m h ± ^ 90 GeV. The same is also true in the constrained m™ ax scenario where /i = 200 GeV 
and A t = — a/6Msusy + fJ>cot/3 are fixed as was already noted in [17] . However, at difference 
to the constrained m™ ax scenario there are regions with larger m ai that are allowed also for 
m h ±< 120 GeV. 

Looking at the [m h ± ) tan f3] -plane one clearly sees the constraint from B u —} r + u T for in- 
termediate tan/5 (show as blue points). For larger tan/5 there is a cancellation between the 
SM and /i ± contributions, which makes this region allowed by B u —> r Jr v T1 but instead the 
constraints from B s —> come into play (red points). It should be noted that the points 

excluded by B u —> t+v t are plotted on top of the constraints from B s —> Similarly 
the constraints from searches for supersymmetric particles are plotted (in yellow) on top of the 
constraints from S-decays, but with the cut \± > 125 GeV there are hardly any points excluded 
by this constraint. Finally we note that for tan/3 < 15 there are allowed points in parameter 
space for essentially any value of m h ± > 90 GeV. 

Turning to the [m ai , tan f3] -plane we see that the region m ai ~ 2ms is allowed by all con- 
straints up to tan/3 < 10, whereas for larger tan/3 there are points allowed by direct Higgs 
boson searches but not allowed by B s fi + ji~ and B u — >> r Jr v T . Given the uncertainties re- 
lated to the indirect constraints from S-decays we conclude that there is a region in parameter 
space with m ai ~ 2ms, m h ± G [90, 170] GeV, and tan /3 G [1, 60] that should be searched for by 
the ATLAS and CMS experiments. It should also be noted that values both above and below 
the threshold m ai = 2m b are allowed by the constraints. 

Before turning to the signal of interest, i.e. Br(h ± —> a{W), we also show in Fig. [l] the 
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FIG. 1. Correlations between the parameters of the scan in the scenario under consideration (see text 
for details) with the various constraints applied as follows. All points in black are viable model points 
but are excluded by HiggsBounds. The coloured points are allowed by HiggsBounds and plotted in 
the following order: red points are excluded by B s —> blue points are excluded by B u —> r + z/ r , 

yellow points are excluded by direct searches for supersymmetric particles, and finally green points 
are allowed by all constraints considered. 

effects of the various constraints when projected onto the [ft, A] and \A t , fj] -planes. From the 
first of these plots one clearly sees the constraint \J k 2 + A 2 < 0.7 which arises from requiring 
perturbativity up to the GUT scale. From the second we see that the constraints imply \A t \ < 
3500 GeV, which essentially follows from the radiative corrections to the lightest CP-even Higgs 
becoming small or even negative for large \A t \ relative to the value Msusy = 1 TeV that we 
are using. We also see that for \i there are hardly no experimental constraints in the region 
considered. On the other hand, if we would extend \i to values smaller than 125 GeV then all 
those points would be excluded by searches for supersymmetric particles. 

As promised we turn now to the resulting branching ratios for the decay /i ± —> a{W '. From 
Fig. [2] we observe the following general feature, the branching ratio can be large as soon as 
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FIG. 2. The branching ratio for h —> a{W when 
consideration (see text for details) with the variou; 
as in Fig. [T] 
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scanning over the parameters in the scenario under 
s constraints applied. The colour coding is the same 



the channel is open (m ai < m^± — mw) except for large tan/3 where the decay /i ± —> rv T 
becomes dominant. Concentrating on those points that pass all the constraints considered and 
the region m ai ~ 2rriB we also see from the lower right plot that Br{h ± —± a\W) > 0.9 as long 
as tan /3 < 10 and m h ± > 100 GeV. Thus we can conclude that not only is the parameter space 
region m ai ~ c 2rriB ) rrih± G [100, 170] GeV, and tan/3 G [1, 10] allowed - in this region the decay 
a\W is also dominant. In the next section we will exemplify how the so far unexplored 
region of parameter space with m ai G [10, 12] GeV can be probed by searching for — >• a\W 
with ai bb in tt-production at the LHC. 

Before ending this section we also show in Fig. [3] the branching rations for the decay chains 
of interest, i.e. t —> t —> bh + —> ba\W^ and t —> bh + —± ba{W —> bbbW as a function of 
tan f3 when restricting to parameter space points with m^± < 160 GeV. As can be seen from the 
figure, points with Br(t —> bh + ) as large as 0.2 are still allowed when the decay /i ± —± a{W is 
included. This should be compared with the experimental constraints from ATLAS and CMS 
which so far have assumed that Br(h ± tv t ) = 1 giving a limit Br(t — )■ bh + ) < few percent 

[23 ESI. 



IV. SEARCH STRATEGY 



In the following section we will perform a signal-to-background analysis for three different 
charged Higgs masses: = 100, 130 and 150 GeV respectively. For definiteness we have 

used tan /3 = 50 when simulating the signal but the end results will not depend on this value. 
The mass of the a\ is set to 11 GeV throughout as an example of a small a\ mass that is just 
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FIG. 3. The branching ratios for the decay chains t — » 6/i + , t bh + — » baiW, and t — > 6/i + — >► 
baiW — >> 666VF respectively when scanning over the parameters in the scenario under consideration 
(see text for details) with the various constraints applied and only considering points with m h ± < 160 
GeV. The colour coding is the same as in Fig. [T] 



above 66-threshold. The three charged Higgs masses are chosen to illustrate different kinematic 
properties: at m h ± = 150 GeV the 6- jet from the t —> bh + decay will be rather soft whereas 
the a\ from the h + will be harder. For = 100 GeV the situation will be the opposite, and 
in the intermediate case rrth± = 130 GeV both jets can be relatively hard and in addition the 
available phase-space will be largest. 

We aim to reconstruct the signal process where one of the t-quarks decays leptonically via 
a{W with W —± tvg and the other hadronically via W —± jj as illustrated in Fig. [4j 
All cross-sections have been corrected for these enforced W decays as well as a factor of 2 for 
taking account of the process also with interchanged roles between the t and the i. Because 
the a\ is supposed to decay close to threshold to 66, we aim for a reconstruction where the two 
6's from the a\ are clustered together to give a single 66-jet. 

As backgrounds to the process we consider the irreducible tibb as well as, because of its 
higher magnitude, ti with one jet being accidentally 6-tagged (weighted with a mis-tagging 
probability, assumed to be 0.01 (29j 130]). In order to include also the single top contributions 
to the background we have simulated the processes ibW + bb and tbW + respectively, but in the 
following we will denote them as tibb and ti for simplicity. For other reducible backgrounds, 
such as W + 6 jets, we assume that similar procedures can be applied as in the ti cross-section 
determination. For example, requiring two 6-tagged jets reduces the + 6jets background to ti 
production to about 10% [31J. Several cuts are applied to strengthen the signal and to suppress 
the background as will be discussed in the following. 

A center-of-mass energy of 8 TeV at the Large Hadron Collider is assumed throughout the 
whole analysis. For the generation of the hard matrix elements, we use MadGraph 5 [32] with a 
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FIG. 4. Illustrations of tt production in proton-proton collisions with the subsequent decay chain 
considered in this paper. 

fixed renormalization and factorization scale (set as default to the Z mass) and the 'CTEQ6LT 
parton distribution functions. To supply MadGraph with the proper parameters of the signal 
we have for simplicity used a simple two Higgs Doublet Model with the masses given above as 
implemented in the Two Higgs Doublet Model Calculator 2HDMC [33J. The only important 
difference to the NMSSM then arises from the a{W decay giving and extra factor cos 2 6 a 

(cf. Table [I]). All other steps to generate complete events, such as radiation, underlying events 
and hadronization, are carried out using Pythia 8. We start from bare samples of 100000 
events for the different signals as well as the tibb background whereas for the tt background we 
have 50 times higher statistics to start with. There is no detector simulation included in this 
exploratory analysis but we have simulated 6-tagging in a simplified way as detailed below. 

For all processes we use the leading order cross-sections obtained from MadGraph. On the 
one hand this means that there is an overall scale factor which is more or less the same for 
both signal and background. On the other hand the rates will be lower than what would have 
been result if higher order cross-sections had been used. All in all this means that the signal 
over background rates we find will be underestimated in this respect. For example we get a LO 
cross-section for pp —> tt process in the SM of 138 pb to be compared with the NNLL resummed 
result of 232 pb [S]. 



A. Reconstruction of the leptonic W 

To reconstruct the leptonic W ) we first need to identify the charged lepton (e or ji) associated 
to the hard process. The transverse momentum (pj_) and pseudo-rapidity (^-distributions 
are shown in Fig. [5| After applying cuts on the lepton kinematics (p± > 20 GeV, \rj\ < 
2.5), we require the summed p± of the surrounding particles in an (77, 0)-cone of size Ai? = 
yj Arj 2 + Ac/) 2 = 0.3 around the lepton to be less than 10 GeV to call it isolated^} On the whole 
event then, we require to have precisely one isolated lepton in the final state. 

The next step in reconstructing the leptonic W is to identify the missing energy (MET) /missing 
p± of the event with the transverse momentum of the neutrino. Assuming the W to be on 
mass-shell and using a mass-less neutrino then leaves two possible solutions for the longitudinal 



The numerical values for the p± cut and the cone size have been optimized by observation of the changes in 
efficiency and purity when varying the cuts 
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FIG. 5. Transverse momentum and pseudo-rapidity of the charged lepton from the W decay and 
missing transverse energy (MET) before applying cuts. 



momentum of the neutrino p vz = A ± sjD with 



A = 



m 2 w vu + 2p lz (p v± pi ± ) 



D = ^L(^f - (p v± f(p lA _) 2 + (p v± p lA _f + (p„ ± p l± )m 2 w ) 



(6) 
(7) 



In order to have a more accurate reconstruction of p uz we have applied a cut f]j_ > 30 GeV. 
In addition, different selection criteria for the choice of the sign have been examined, e.g. the 
invariant mass of the reconstructed h ± . Among these, the most viable one turned out to be 
a simple selection of the smaller \p uz \, which is correct in roughly three quarters of the signal 
events. 



B. Jet reconstruction and tagging 

For the reconstruction of the hadronic part of the event, we consider different jet clustering 
schemes (anti-kT [35], Cambridge/ Aachen [36], kT j37J|38]) as well as cone-sizes. All particles, 
except neutrinos and the isolated lepton, in the rapidity range \rj\ < 4.9 are fed into Fast Jet [39J 
[40] . The resulting clustered jets are required to have p± > 20 GeV. Afterwards, a simplified 
b-tagging is simulated for all jets in the region \rj\ < 2.5 by comparing the (77, 0) of the jets to 
the 6-quarks of the hard process. All jets within AR = 0.4 are then classified as 6-jets. 
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W reconstruction Top reconstruction 




FIG. 7. Mass reconstruction of the hadronically decaying W (left) and t-quark (right) for the signal 
with m h ± = 130 GeV when R = 0.5 (left) and using the anti-kT algorithm (right). 



The jet-algorithms requires to specify the distance measure R used when calculating the jet 
measure. Since we want to cluster the bb pair from the a\ into one jet, it is crucial that the 
distance between the two subjets is not too large. Fig. [6] shows the distance measure between 
the two 6-quarks on part on level for different a\ masses. As can be seen from the figure, for 
reasonable clustering cone sizes (around 0.4 to 0.6) and with m ai in the region of interest for 
this analysis, the two 6-quarks will most likely be clustered together as a single jet. 

As an ideal reconstruction will now give rise to three 6-jets, the correct reconstruction of the 
/i ± will be enhanced by the identification of the "wrong" b-jet which comes from the t together 
with the hadronically decaying W~ . The strategy to achieve this here is to first find the pair of 
untagged jets that is closest to the W mass and then combine this pair with the 6-jet that gives 
a mass closest to the t mass. This 6-jet will then be excluded in the /i ± reconstruction, reducing 
the number of 6- jets which have to be considered (in an event with a so far correct clustering 
in the desired way) to two. In addition to this, we put cuts on the quality of the reconstruction 
by requiring the reconstructed masses shown in Fig. [7] to be in the regions mw ± 20 GeV or 
m t ± 30 GeV respectively. We have checked that the reconstruction of the W and t masses are 
quite independent of the choice of jet clustering scheme and cone-size as is also shown in Fig. [7| 
Only the 0.3 cone size gives a slightly inferior top reconstruction, but then the a\ will also not 
give rise to a single 6-jet. If not mentioned differently, we thus use the anti-kT algorithm with 
R = 0.5 in the following. 
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FIG. 8. Jet (left) and 6- jet (right) multiplicity for the signal and backgrounds when using the anti-kT 
algorithm with R = 0.5. 
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FIG. 9. Reconstructed charged Higgs masses, using either the higher-p^ (left) or the lower-p^ (right) 
6-jet combined with the leptonically decaying W. 



Using the jet-reconstruction outlined above we obtain the jet and 6-jet multiplicity respec- 
tively shown in Fig. [8| As is clear from the figure, the number of jets peaks around the expected 
five and the number of 6-jets is a typically smaller than ideal which is due to the limited 6- 
tagging region enforced. Also note the large decrease in 6-jet multiplicity for the tt sample from 
2 to 3. Here the 3 6-jet sample arises from gluon splitting into 66 and we do not take it into 
account below since that would amount to double counting with the tibb background. 

Given the large background from tt for the two 6-jet sample we resort to requiring at least 
three 6-jets. Assuming that one of the 6-jets has been identified as coming from the t-quark 
this leaves two 6-jets that may originate from the h ± . Due to the strong dependence of the 
hardness of the 6-jet from the t — >■ h + b decay on the mass, we consider both of these 
remaining solutions in general as possible candidates for the reconstruction. The resulting 
distributions when combining with the leptonically decaying W are shown in Fig. [9j In the 
/i ± reconstruction, we thus require 3 b-jets for the signals and the tibb background, while we 
require 2 6-jets in the tt sample. For the latter we then assume that any of the non 6-jets inside 
the b-tagger region (\r\\ < 2.5) can be mis-tagged with a probability of 0.01 per jet |29l [30] . 
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C. Signal significance and reach 



To estimate the signal reach as well as its significance, we choose a common window of 90 to 
160 GeV to integrate the signal as well as the backgrounds. We correct for 6-tagging efficiencies 
which we assume to be 0.6 per 6-jet. The resulting cross-sections are shown in Table [TTJ 





Vpeak [fb] 


signal /background 


100 GeV 


130 GeV 


150 GeV 


tibb 


tt 


£ BG 


High p± channel 


0.78 


7.9 


3.7 


3.0 


0.9 


3.9 


Low p± channel 


0.97 


9.2 


4.7 


4.8 


1.2 


6.0 



TABLE II. Integrated signal cross-sections (cf figure [9]) for the three different charged Higgs masses 
as well as the two backgrounds. The integration region is 90 to 160 GeV 



S/\/B ratios for the three different mass cases are then calculated for an integrated lumi- 
nosity of 20 fb -1 and summarized in Table III. The table also shows the branching ratios at the 
simulated parameter points, as well as the extrapolations to branching ratios necessary for a 5a 
discovery. From the table it is clear that for m h ± ~ 130 GeV, the discovery reach is maximal. 
For smaller or larger m^±, the limitations in phase space will reduce the branching ratio for 

—> a{W and t h^b respectively. In the former case this means that the standard decay 
channel /i ± rv T will also be significant. 





100 GeV 


130 GeV 


150 GeV 


6-jet selection 


high p± 


low p± 


high p± 


low p± 


high p± 


low p± 


s/Vb 


1.77 


1.76 


18.0 


16.7 


8.5 


8.5 


BR(t -»■ bh + -»■ baiW -> bbbW) 


0.0051 


0.022 


0.015 




0.014 


0.014 


0.0060 


0.0065 


0.0085 


0.0085 



TABLE III. The S/y/B ratios obtained for an integrated luminosity of 20 fb -1 for the different signals 
considered together with the branching ratios for the total decay chain t — >► bh + — » ba{W bbbW 
of the respective parameter points and a linear extrapolation to the critical branching ratio necessary 
for a 5a discovery. 

Before ending this section we note that similarly to the standard decay modes of the charged 
Higgs boson it should be possible to use the spin-correlations between the decay products of 
the two top quarks as a way of enhancing the signal [4TT443] . 



V. COMPATIBILITY WITH POSSIBLE HIGGS SIGNAL 



Recently the ATLAS and CMS experiments announced the combined results of the SM 
Higgs (0) searches 5 fb 1 of data from 2011 [HI [45]. In short they can be summarized 

as follows: the CMS experiment has ruled out the region m$ G [129, 600] GeV at the 95 % 
confidence level as expected whereas in the region [118,129] GeV they have not been able to 
make any exclusion at all even though they expected to be able to do so at the 95 % confidence 
level, the ATLAS experiment has similarly ruled out the regions G [112.9, 115.5] GeV, 
G [131,238] GeV, and G [251,466] GeV at the 95 % confidence level whereas they had 
expected to rule out the region [125, 519] GeV. 
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Instead, both experiments have found an excess of events in the regions ^126 GeV 
and 171^ ^124 GeV for ATLAS and CMS respectively, which when corrected for the so called 
look-elsewhere-effect and combining the different channels has a statistical significance of about 
2a for each of the two experiments. Although not statistically significant these results have 
stirred a lot of excitement in the high-energy physics community j2U H6H5^] . It will most likely 
not be possible to draw any final conclusions about whether this is a true signal or not before 
the end of the run in 2012 which will give another ~ 20 fb _1 of data. 

One of the most important properties of the possible signal at the LHC is that the (j) —> 77 
channel is similar to what is expected from the SM. Therefore we start by considering the would 
be signal from the h\ as well as the h 2 compared to what is expected in the SM for a Higgs 
boson {(j)) with the same mass. Assuming that gluon-gluon fusion dominates the production, 
which we have verified is always the case in the scenarios we consider, this ratio is given by 

Rhi a (99 Mnmssm Br (hj ->> 77)nmssm ^ T(hj -)> ^)nmssm Br (hi -> 77)nmssm 
9911 ~ o(gg 0) SM Br((f) 77) SM " T(</> gg)sM Br ((f) 77)sm ^ ) 

where in the second equality, following for example [5H [53], we have made the implicit as- 
sumption that the difference in radiative corrections for the production and decay processes are 
canceled in the ratio. 



The results obtained for this quantity when using the same scan as in section III are displayed 



in Fig. 10. From the figure it is clear that if the possible signal seen at the LHC is the hi then 
we would have to have m ai > 60 GeV. Even so Rgg 71 does not become larger than ~ 0.5 
for m/ ll > 120 GeV in our scan and there appears to be an upper bound m/ ll < 122 GeV. 
However, the difference in mass compared to the possible observation is so small that it should 
be considered to be within the theoretical uncertainty [^J In any case it is clear that it is hardly 
possible to have a light ai if it is the hi that has been seen at the LHC. 

Next we turn to the possibility that it is the h 2 that has been seen by the LHC experiments. 



The results of the scan are also shown in Fig. 10, As can be seen from the figure the results 
are more promising in this case. There are points with i?^| 77 ^0.1 also for small m ai and 
intermediate m^± that have rrih 2 in the region indicated by the LHC experiments. 



In order to explore this possibility more closely we show in Fig. [TT] the results when fixing 
the ai mass to 11 GeV as was done in section [TV] for the three cases m h ± = 100, 130, 160 GeV. 
As can be seen from the figure, for the intermediate charged Higgs mass it is possible to reach 
Rggjj ~ 0.15. However it should be noted that in this case, as also shown in the figure, the 
branching fraction for a{W is quite small meaning that the standard decay channel 

— >> rv T can be used even though with a slightly reduced branching fraction. 

The last logical possibility would be that it is the /i 3 that has been observed by the CERN 
experiments. However, in the scenarios we consider it turns out that R^ gil is always small. 



The trends seen in Fig. 11 can be understood on more general grounds from the difficulties 
of having a light a\ with mass m ai < and at the same time have a large Rg z gil - The 

first problem is that unless the triple Higgs coupling gh iaiai is small the decay hi —> a\a\ will 
become dominant. Looking at the structures of gh iai an which for example can be found in [8], 
this means that both A and n typically have to be small. Secondly the decay hi — >■ a\Z will 
also dominate if m ai < — mz unless the corresponding reduced coupling given in table [I] 
is small. In other words we need (cos/3S^i — sin/3S^ 2 ) cos #,4 to be small. There are essentially 
three ways to achieve this. If cos #,4 is small then a\ will be mainly singlet like and decouple. 
However, then a\W will also be small. The second possibility is that S^i and S^ 2 are 



2 For example, including the full one- loop corrections and the two-loop ones from the t and b Yukawa couplings 
could push the hi mass bit higher in compliance with the possible signal. Similarly increasing Msusy would 
also increase • 
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FIG. 10. The signal for hi 77 (left) and /12 — > 77 (right) relative to the standard model when 
scanning over the parameters in the scenario under consideration (see text for details) with the various 
constraints applied. The colour coding is the same as in Fig. [T] 



small, but then the would-be-signal would be mainly singlet-like and not produced in the first 
place, so S^ 3 has to be small. Finally then, the combination (cos f3Sn — sin /3S^ 2 ) could be small 
but then the complementary combination (sin/3S^i + cos/5S^ 2 ) would have to be large giving 
an increased coupling hi VV. All in all this means that is difficult to have a light a\ that is 
not decoupled and still have a large ife™, although it cannot be completely ruled out. 



VI. SUMMARY AND CONCLUSIONS 



We have considered a well motivated class of supersymmetric extensions of the standard 
model with non-minimal Higgs sector - namely the CP-conserving NMSSM. In these types 
of models the additional Higgs singlet can modify the phenomenology of the Higgs sector in 
many different ways. In this paper we have specifically addressed the possibility of having a 
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FIG. 11. The signal for /12 — > 77 relative to the standard model (left) and the branching fraction for 
a{W (right) in a scan with m ai — 11 GeV and m h ± fixed at 100 (top), 130 (middle), and 160 
(bottom) GeV respectively. The colour coding for the various constraints is the same as in Fig. [I] 



light CP-odd Higgs boson close to, but still above the bb threshold. This in turn means that 
the light charged Higgs boson, with mass m h ± < m t , can decay into a\W in addition to the 
standard decays — >■ rv T) thus invalidating the interpretations made of charged Higgs bosons 
searches assuming that Br(h ± — >> tv t ) = 1. 

When investigating the viability of these types of scenarios we have found that the experi- 
mental constraints from direct searches are quite weak even when taking the latest constraints 
from LHC into account. The constraints from indirect searches in £>-decays are more constrain- 
ing but also more model dependent. Even when including the results from the most important 
ones, namely B u — >> rv T and B s —> we still find a region of parameter space that is 

allowed with tan/3 G [1,10]. This is precisely the same region as the one where the decay 

—> a{W can be dominant. 

The phenomenology of these types of scenarios is special in that, due to its low mass, the 
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di will decay into a single bb jet. Even so we have shown that it is possible to reconstruct the 
decay —> a\W using standard jet finding algorithms when the W decays leptonically. This 
requires to use the missing transverse momentum to calculate the four momentum of the W, 
which can then be combined with the ai-jet to give a mass peak at m h ±. The other t quark 
is assumed to decay hadronically according to t —> bjj giving an additional handle to identify 
the events of interest. The most important background is thus the irreducible one from tibb 
production but we have also taken into account the tt background by considering the possibility 
of mis-tagging ordinary jets as 6-jets. 

Based on our study we find that with an integrated luminosity of 20 fb _1 it should be 
possible to discover a charged Higgs bosons in these types of scenarios as long as the combined 
branching fraction for the decay chain t — >> bh + — »> ba\W —± bbbW is larger than « 0.01. 

Finally we have also investigated the phenomenological consequences of the possible Higgs 
signal seen at the LHC on the types of scenarios we consider. We find that it is difficult to have 
a light a\ that is not decoupled and at the same time have a combined production and decay 
into 77 for one of the CP-even Higgs bosons with mass ^ 125 GeV. As a consequence we have 
not been able to find regions of parameter space where the — >> a\W decay dominates and 
at the same time are compatible with the possible Higgs signal. This means that irrespectively 
of whether the possible Higgs signal is substantiated or not, the LHC experiments should be 
able to either discover or put very tight constraints on a light charged Higgs bosons also in the 
NMSSM. 

Note added: 

On July 4, 2012 the ATLAS and CMS experiments announced the discovery of a new Higgs-like 
particle in the same mass region as the previous observations already cited in the text. 
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